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Abstract. Reflection of an obliquely incident solitary wave onto a vertical wall 
is studied analytically and experimentally. We use the Kadomtsev-Petviashivili 
(KP) equation to analyze the evolution and its asymptotic state. Laboratory ex- 
periments are performed using the laser induced fluorescent (LIF) technique, and 
detailed features and ampliflcations at the wall are measured. We flnd that proper 
physical interpretation must be made for the KP solutions when the experimental 
results are compared with the theory under the assumptions of quasi-two dimen- 
sionality and weak nonlinearity. Due to the lack of physical interpretation of the 
theory, the numerical results were previously thought not in good agreement with 
the theory. With proper treatment, the KP theory provides an excellent model 
to predict the present laboratory results as well as the previous numerical re- 
sults. The theory also indicates that the present laboratory apparatus must be 
too short to achieve the asymptotic state. The laboratory and numerical results 
suggest that the maximum of the predicted four-fold ampliflcation would be diffi- 
cult to be realized in the real-fluid environment. The reality of this amplification 
remains obscure. 



1 Introduction 

Due to the nonlinearity effect of surface waves, the interference of multiple waves causes local 
amplification. In deep water, anomalously large waves are often called freak (or rogue) waves. 
The cause of this phenomenon is not fully understood, but has been speculated by various 
theories. A traditional explanation is the role of currents. When strong currents meet waves 
moving in the opposite direction, the currents focus and concentrate sets of waves, causing 
individual peaks that are exceptionally higher than their surroundings. The amplification can 
also occur when disparate trains of waves meet together. Such intersections sometimes cause 
waves to be higher and steeper. 

Freak-wave formation in shallow water was analyzed by Pelinovsky, et al. |24| within the 
framework of nonlinear-dispersive theory in single propagation space. They attributed the main 
mechanism of freak-wave formation to the spatial-temporal focusing of radiative wave packets 
and a solitary wave due to the difference in their propagation speeds. They pointed out that a 
freak wave is rare, short-lived, and always has a small ratio of nonlinearity to dispersion. 

Wave refraction, reflection and diffraction lead to obliquely interacting waves in shallow 
water, and when their superposed amplitudes are sufflciently large, nonlinear effects can have 
striking effects on the resulting surface wave patterns. The nonlinear interaction of obliquely 
propagating solitary waves in similar directions can cause a high wave hump resulting from 
the crossing of the waves in the two-dimensional propagation space. Wiegel f3T] described that 
regular reflection of a solitary wave at a vertical rigid wall does not occur for sufficiently small 
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Fig. 1: Definition sketch for Macli reflection: ipo, incident wave angle. 



angle of incidence, being replaced by a Mach reflection (a geometrically similar reflection from 
acoustics ^2J). The apex of the incident and reflected wave separates away from the wall and 
is joined to it by a third solitary wave that perpendicularly intersects the wall as depicted 
in Figure [l] Wiegel [3T] further reported that when the 1946 Aleutian Tsunami hit the town 
of Hilo, Hawaii, a Mach-stem phenomenon was observed along the cliffs forming the western 
boundary of the Hilo Bay (Figure [2]) . 

Here we focus on such shallow-water wave interactions: two identical solitary waves propa- 
gating with a small oblique angle 2-00 • This condition is equivalent to the Mach reflection with 
the incident wave angle tpQ to a perfectly reflective vertical wall. Miles [27, presented a the- 
oretical analysis of the Mach reflection at its asymptotic state. His theory predicts the extraor- 
dinary four-fold amplification of the stem wave. There were several laboratory and numerical 
experiments to attempt to validate Miles's four-fold amplification, but with no definitive success 
(for numerical experiments, e.g.l8],[27],[T5],[T], and for laboratory experiments, [25], [20]). 

In this paper, we reanalyze Mach reflection in a systematic fashion based on the Kadomtsev- 
Petviashvili (KP) equation. Our theoretical analysis allows us to predict not only the asymptotic 
states but also the developing states of the Mach reflection phenomenon. Then, we present the 
laboratory data obtained from the precision experiments, and quantitatively discuss agreement 
and discrepancies between the theoretical predictions and the laboratory results: also reanalyze 
the previous numerical results by Tanaka 27J. 



2 The KP solitons 



Based on the weakly nonlinear, weakly dispersive waves having the small oblique angles ipQ = 
C(e) with < e <C 1, we assume the following orders in the wave amplitude oq and the constant 
water depth ho at the quiescent state: 



"0 X / ho 



2 



0{e), (^gj =0(e), tan2V'o = 0(e), 

where Aq is the wavelength scale. Then, the Euler formulation results in the Kadomtsev- 
Petviashvili (KP) equation of the form 

, 3co cohl ^ , Co „ 

Vt + coVx + + —^-n^xx J + -jVyy = 0, (1) 
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Fig. 2: Wave in the form of a bore, approaching Wailuku River, Hilo, during the 1946 tsunami. 
It appears to be moving along the seaclifF as a Mach stem. (Courtesy Modern Camera Center, 
Hilo). 



in which rj is the water-surface elevation from ho, and cq 
equation for one-line soliton can be written as: 



/gho. An exact solution of the KP 



rj = apsech 




. CLq 1 2 / \ 

x + yt&ntpo - Co [ 1 + ^ + - tan ^Jq ] t 



(2) 



Except for the case ^po — 0, the KP solution ^ does not satisfy the KdV equation; hence it does 
not support the KdV soliton. One should also note that the solution ^ becomes unphysical 
for a large angle. This is due to the assumption of the quasi-two dimensionality of the problem. 
The effect of the finite angle can be imposed in the following way: first write the KP solution 
^ in the form, 



rj — flgsech 



4/i3 cos^ i/^o 



X cos "00 + y sin ^0 — Cg cos ■(/'o ( 1 + TTT^ + tan^ -00 ) ^ 

2/iQ 2 



(3) 



Define x '■— a;cos'!/'o + ysin-^o (the normal direction to the one line-soliton wave crest), and 
note cos V'o = + tan^ t/jq = 1^5 tan^ -00 + C(e^)7 we have 



r] = Oosech 



3ao 



4/ig cos^ ipo 



(4) 



Note here that the width of the soliton depends on the angle in the higher order, which is 
evidently unphysical. So we now define an equivalent r] including a higher order correction to 
the quasi-two dimensionality. 



fj — flosech^ 




X - Co 1 



ao_ 
2ho 



with flo ~ 



ao 



cos^ ipo 



(5) 



which gives a better approximation than the KP solution ([2| for the solitary wave with the 
propagation angle i/jq in physical situation. Notice that with this new amplitude ao, the soliton 
solution ([5]) satisfies the KdV equation in the form up to 0{e^), 



3co 



cohl 



Vt + cofjx + T^VVx + -^Vxxx = 



2ho 



(6) 
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The point here is that the amphtude of the soliton observed in the experiment (or numerical 
simulations of the Euler equation) should be represented by gq instead of oq of ([2|. 

As for the KP limit, the solutions expressed by (2]) and ^ are equivalent. None the less, 
it makes a significant difference when the theoretical predictions are compared with the data 
from laboratory and/or numerical experiments — some of the experiments, including the results 
presented in this paper were performed with large amplitude ao/ho and large incident angle 



2.1 Exact soliton solutions of the KP equation 

Let us write the KP equation ([I]) in the following standard form, 

{AuT + 6uux +uxxx)x + ^'^YY ^0, (7) 
where the non-dimensional variables {X, Y, T) and u are defined by 

V = ^ u, X- Cot = hQX, y = haY, ^ T. (8) 

6 2cq 

Notice that the one-line soliton solution in the standard form is given by 

u = Aosech^ {X + Y tan V'o - CT) , (9) 

where C = ^Ag + | tan^ Vo and the amplitude Ao has the relation with the physical variables, 

3ao 3ao 2 / 
^" = 2^ = 2^ '^"^ 

With the standard form of the KP equation in Q, we consider the solution in the form 

u{X,Y,T) = 2[lnT{X,Y,T)]xx- (H) 
where t is referred to as the tau function defined in the Wronskian determinant, 

T = Wr(A,/2,...,/Ar). (12) 

The functions /i's satisfy the linear equations, fy — fxx, fr = —fxxx (see e.g. Il2j). For our 
soliton solutions, we choose the solution form of: 

M 

/z = E a^Je'>^ , 0, = KjX + K^r - K^^T. (13) 

Here the coefficient matrix A = {uij) is a constant N x M matrix, and Kj are constants with 
the ordering ki < • • • < km- Thus each solution u(X,Y,T) is then completely determined by 
the A-matrix and the K-parameters. 

Single soliton solution is obtained in the case with = 1 and AI = 2: We have r — e^' +e^^ 
for Ki < Kj, and u = 2{hiT)xx gives 

u^^{K,-K,)hech'^{0,-0,). (14) 

The line 9i — 6j represents the ridge of the soliton, and we refer to this Zme-soliton as 
soliton with i < j (i.e. Ki < Kj). The amplitude j] and the inclination tan-f/jji jj of this soliton 
are given by 

Ai,i] = ~ '**)^' tanV'iij] = Ki + Kj, (15) 

where — 7r/2 < < 7r/2 is the angle of the line soliton measured counterclockwise from the 
y-axis. This angle also represents the propagation direction of the line soliton (see Figurdsl). 



Will be inserted by the editor 



5 



2.2 The Mach reflection and the KP solutions 



For the Mach reflection problem, we take N — 2, M = 4. Using the Binet-Cauchy theorem for 
the determinant, the r-function r = Wr(/i, /2) can be expressed in the form, 

E ahj)Eii,j), (16) 

l<i<j<4 

where is the 2x2 minor consisting of i-th and j-th columns of the A-matrix, and 

E{hj) — Wr(e^' , ) = {nj — Ki)e^^^^^ . For the regular solutions, we require all of these minors 
to be non-negative (note E(i,j) > with the order Ki < Kj). 
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Fig. 3: Contour plots and the chord diagrams of soliton solutions. Left: 0-type solution. Right: 
(3142)-type solution. Each denotes the [z, j]-soliton. The length of [1, 4]-soliton changes in 
t. The upper (lower) chords represent the asymptotic solitons in y ^ {Y ^ 0). The thicker 
chords correspond to the solitons in the right side {X 3> 0), and the arrows on the chords show 
the pairings in the permutations. 



As was shown in each r-function (16) generates a soliton solution which consists of 

at most two line-solitons for both y ±oo. We consider the following two types which are 
relevant to the solutions of the initial value problems for the Mach reflection: one consists of 
two line-solitons of [1, 2] and [3, 4] for both |y| ^ 0, which is called 0-type soliton ("O" stands 
for original, see [T5],2|); the other one consists of [1,3] and [3,4] line-solitons for Y ^ and 
[1,2] and [2,4] line-solitons for Y <^0. Let us call this sohton (3142)-type, because those four 
line-solitons represent a permutation n = {}^^^*) . It was shown in general (see [3],|4]) that 
each exact solution generated by the Wronskian determinant with ( |13[ ) can be parametrized 
by a unique element of the permutation group. In this representation, 0-type is expressed as 
(2143)-type solution. Figure 3] illustrates the contour plots of 0-type and (3142)-type solutions 
in the XY-p\a,ne, and the corresponding chord diagrams which represent each soliton as a chord 
joining a pair of k^'s following its permutation representation. The upper chords represent the 
asymptotic solitons [i,j] for Y ^ 0, and the lower chords for the asymptotic solitons [i,j] for 
y <C 0. Note also that the length and the location of each chord give the amplitude and the 
angle of the corresponding soliton (cf. ( 15 1). The A-matrices for those solutions are respectively 
given by 

^o^(o01&)' ^(3142) = (Jo5 ^"^l , (17) 
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where a,b,c > are constants determining the locations of the soUtons (see [4]). Notice that 
the r-function for the (3142)-type contains five exponential terms, 



r = (k3 - Ki ) 4_ (k4 - Ki ) 6 e^i + (k3 - K2 ) a e 

+(K4-K2)a6e*^+^^ + (K4-K3)ce^^+^% 



(18) 



and the r-function for 0-type with c = in ( 18 ) contains only four terms 



Let us fix the amplitudes and the angles of the solitons in the positive X regions for both O- 
and (3142)-types, so that those solutions are symmetric with respect to the X-axis (see Figure 
©: 

-4[i.2] = ^[3,4] (f™ O - type) 
^[1,3] = ^[2,4] (for (3142) - type) 

(19) 

-V^[i,2] = V'[3,4] > (for O - type) 
-V'[i,3] = i'[2A] > (for (3142) - type) 



An 



Then from (151, one can find the K-parameters in terms of Aq and tan ■00 with ki 
K2 = — K3 (due to the symmetry): In the case of 0-type, we have 



Ki = - i (tanV'o + 



K2 = (tan -00 



-K4 and 



(20) 



The ordering K2 < K3 then implies tan-0o > -\/2^o- On the other hand, for the (3142)-type, we 
have 

Ki = - - (^tanV'o + V2A)) , ^'^=2 {^^'^'^^ ~ ■ (21) 

The ordering K2 < K3 implies tan^o < V^Aq. 

Thus, if all the solitons in the positive X-region have the same amplitude Aq for both O- 
and (3142)-types, then an 0-type solution arises when tan -00 > V^Aq, and a (3142)-type when 
tan -00 < V^Aq- Then the limiting value at K2 = K3 {— 0) defines the critical angle ipc, 



(22) 



Note from (18) that at the critical angle, i.e. K2 = K3, the r-function has only three exponential 



terms, and this gives a "Y" -shape resonant solution as Miles noted j22j . 

One can also see from (15) that for (3142)-type solution, the solitons in the negative X- 
region are smaller than those in the positive region, i.e. ^[3,4] = ^[1,2] = 5 tan^ ipo < Aq, and 
the angles of those in the negative x-regions do not depend on ipg and "0(3, 4] — ~^[i.2] = "00- 
Two sets of three solitons {[1, 3], [1, 4], [3, 4]} and {[2, 4], [1, 4], [1, 2]} are both in the soliton 
resonant state. These properties of the (3142)-type solution are the same as those of Miles's 
asymptotic solution for the Mach reflection of a shallow water soliton [55] . However one should 
note that the [1, 4]-soliton corresponding to the Mach stem becomes a soliton solution only 
in an asymptotic sense. Then the exact solution of (3142)-type can provide an estimate of a 
propagation distance, at which the amplitude is sufficiently developed. 



3 The Mach reflection in shallow water waves 

Miles [5T],[22] considered an oblique interaction of two linc-solitons using O-type solutions. He 
observed that resonance occurs at the critical angle ^pc, and when the angle 0o is smaller than ipc, 
the 0-type solution becomes singular (recall that at the critical angle ipc, one of the exponential 
term in the r-function vanishes) . He also noticed a similarity between this resonant interaction 
and the Mach reflection. This may be illustrated by the left figure of Figure[4j where an incidence 
wave shown by the vertical line is propagating to the right, and it hits a rigid wall with the angle 
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Fig. 4: The Mach reflection. The left figure iUustrates a semi-infinite line-soliton propagating 
parallel to the wall. The right figure is an equivalent system to the left one when we ignore the 
viscous effect on the wall. The resulting wave pattern is of (3142)-soliton solution. The angle (j) 
becomes zero if the initial angle satisfies tpQ > i.e. no stem. 



—ipi^: the incidence wave angle ip^ is measured counterclockwise from the axis perpendicular to 
the wall (see also ^SJ). If the angle of the incidence wave (equivalently the inclination angle of 
the wall) is large, the reflected wave behind the incidence wave has the same angle V'o, i-e. a 
regular reflection occurs. However, if the angle is small, then an intermediate wave called the 
Mach stem appears as shown in Figure |4j The critical angle for the Mach reflection is given by 
the angle ipc- The Mach stem, the incident wave and the reflected wave interact resonantly, and 
those three waves form a resonant triplet. The right one in Figure |4] illustrates an equivalent 
system of the wave propagation in the left figure, (ignoring the effect of viscosity on the wall, 
i.e. no boundary layer). At the point O, the initial wave has V-shape with the angle i/jq, which 
forms the initial data for the simulation. Then the numerical simulation describes the refiection 
of line-soliton with an inclined wall, and these results explain well the appearance of the Mach 
reflection in terms of the exact soliton solution of (3142)-type (see [T7].|14j'). 

The maximum amplitude for this problem occurs at the wall and this can be obtained by 
the following formula (see [12],[1]): For the 0-type solution with tan -00 > tan^c = \/2Ao, 

=, with fc:=^>l, (23) 



1 + Vl - k-^ 
and for (3142)-type solution when A; < 1, 

Ai,4] - ^(«4-Ki)2 = Ao(l + fc)^ (24) 

Note here that at the critical angle V'o = "^c (i-e. k — 1), the amplitude takes the maximum 

"max = ^[1,4] = 4Ao (1.6. = MmaxMo = ^[l,4]Mo = 4). 



4 Laboratory Experiments 

Laboratory experiments were performed in a wave tank (7.3 m long, 3.6 m wide, and 0.30 m 
deep). The wave basin is equipped with a 16-axis directional-wave maker system along the 3.6- 
m long sidewall, capable of generating arbitrary-shaped, multi-directional waves. The paddles 
are made of PVDF (Polyvinylidene fluoride) plates that are moved horizontally in piston-like 
motions. Each paddle has a maximum horizontal stoke of 55cm more than adequate to generate 
very long waves with a water depth of 6.0 cm. While the wavemaker is capable of generating 
arbitrary-shaped, multi-directional waves, we chose to generate a solitary wave in the normal 
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direction (along the sidewalls). An obliquely incident solitary wave was created by placing a 
2.54 cm thick Plexiglas plate vertically at a prescribed azimuth angle from the tanks sidewall 
(see Figure [5]). Because the wave paddles are driven synchronously along the sidewalls, any 
ambiguity associated with potential deviation caused by the paddle deformation is eliminated. 
Solitary waves are generated using the algorithm based on KdV solitons pDj . 



Wave paddles 





Fig. 5: Schematic drawings of the laboratory apparatus: (a) a plan view showing the waveguide 
that creates an oblique wave reflection and the tank frame that can be seen through the bed 
that is made of glass plates; (b) an elevation view. 



To examine temporal and spatial variations of water-surface profiles, we used the Laser 
Induced Fluorescent (LIF) method. Figure |6] shows a setup for the LIF method used in this 
study. A laser beam (a 5W diode-pumped solid-state laser head mounted on the traversable 
carriage) is converted to a thin laser sheet using a cylindrical lens. Two front surface mirrors 
direct the laser sheet to illuminate the vertical plane perpendicular to the wall. With the aid 
of fluorescein dye dissolved in water, the vertical laser-sheet illumination from above induces 
the dyed water to fluoresce and identifles the water-surface profile directly and non-intrusively. 
The illuminated images are recorded by a high-speed high- resolution video camera (1280 x 
1024 pixels and 100 frames per sec). The captured images are rectified and processed with the 
calibrated image so that the resulting images can be analyzed quantitatively. Consequently, 
the LIF method is a highly accurate measurement technique to capture the spatiotemporal 
water-surface profiles in the non-intrusive manner. 

The origin of coordinates is taken at the edge of the oblique vertical wall; x-direction points 
horizontally along the wall; y-direction points perpendicularly away from the wall; and z- 
direction points upward. The present experiments focus on the conditions of the quiescent 
water depth ho = 6.0 cm and the incident wave angles ipo = 20°, 30° and 40°. Hereinafter, we 
present the data in the nondimensionalized forms consistent with (|8| unless stated otherwise. 
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Fig. 6: Setup for the Laser-Induced Fluorescent (LIF) method, comprised of the 5W laser, the 
cylindrical lens, and two front-face mirrors. The water dyed with fluorescein (green) fluoresces 
when excited by the laser sheet. 



CompiUng three LIF segments of the water-surface profiles in the y-direction yields a water- 
surface map in the FT-plane with adequate coverage, < y < 13.0, as shown in Figure]?^. The 
figure represents the wave profile at X = 71.1 for the case of the incident wave amplitude = 
0.212 with the oblique wall V'o = 30°. The plot in Figure[7^ is made from 250 slices of the spatial 
profiles — 100 slices per second — with approximately 3000-pixel resolution in the K-direction. 
As expected, the stem-wave formation is realized, in which the apices of the incident and 
reflected waves separate away from the wall by the third wave that perpendicularly intersects 
the wall. Figure^ shows the corresponding (3142)-type exact solution of the KP equation. For 



Ao = 0.212 _and V-o = 30°, we found from ([T5| that = (-0.614, -0.037, 0.037, 0.614), and the 



A-matrix (17) is given by 

'1 8.797 -1.128 



1 0.530 

This choice of the A-matrix places the incidence wave crossing at the origin {X = 0,Y = 0) 
at T = [4] . The resulting wave profile shown in Figure [TjD is in excellent agreement with the 
laboratory observation, although a careful observation reveals that the measured stem-wave 
amplitude is slightly lower and the reflected wave amplitude is also lower than the prediction. 
This is because the wave reflection in the laboratory is still in the process of being established, 
while the reflected wave is assigned a priori as the [3, 4]-soliton in the KP theory. 

Growing stem- wave amplification factors a {— A[i ,^^/Aq) induced by a variety of incident 
waves with amplitudes 0.086 < Aq < 0.413 and ipo = 30° are presented in FigurejS] No complete 
data are presented for the case of Aq ~ 0.413 because of the wave breaking after X = 50.8. 

The stem amplification continues to grow, and it is evident that the limited physical di- 
mension of the laboratory apparatus prevents the Mach stem from reaching its fully developed 
state. In Table [ij the experimental results of the amplification factor a are compared with 
those KP exact solutions (i.e. 0-type for A; > 1 and (3142)-type for k < 1). At X — 71.1, the 
farthest measuring location in our experiments, the experimental results are in good agreement 
with those calculated from the KP exact solutions. We estimate the asymptotic amplifications 
using the exponential curve fitting to the data: a — ae^^^ + c , where a, b and c are positive 
constants to be determined. Those estimates also agree well, except for the case near k = 1 
(i.e. ^0 = 0.161 and 0.212). Near k = 1, the discrepancy is significant and the theory predicts 
approximately 45% greater than the amplification projected from the experimental results. 
This suggests that while the growth rate be reduced gradually, the theoretical growth must be 
sustained for a much longer time than the actual growth in the real (viscous) fiuid environment. 

The present laboratory data are also compared with the data from Tanaka's numerical 
experiments |27j : his experiments are based on the high-order spectral method. This higher- 
order model allows him to study conditions less restricted in the wave amplitude Aq than the 
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(b) 

Fig. 7: Two views of the temporal variation of the water-surface profile in the y-direction 
(perpendicular to the wall) at X = 71.1: (a) the experimental result, and (b) the corresponding 
(3142)-type exact soliton solution of the KP equation. The incident wave amplitude Aq = 0.212, 
and the angle ipo = 30°. The data were obtained by the LIF method with the vertical laser 
sheet along the F-direction. 



KP limit. Tanaka experienced that it took long propagation to achieve the saturated conditions 
in stem amplitude {X ~ 100 — 300). Like other laboratory or numerical experiments (e.g. 
[H] 1 [H] 1 [1] I [2l] , [20] ) 1 his model fails to simulate the four- fold amplification of the stem wave. 
The maximum stem-wave amplification computed by Tanaka [27, is a = 2.90 aA, k — 1.03 in 
the case of Aq = 0.28. 

Note that Tanaka's maximum amplification a ~ 2.90 is presented at fc = 0.695 instead of 
1.03. It is because Tanaka followed Miles's analysis [22 j which assumes the incident wave angle 
ihn to be small, replacing tan -00 — ^ V'o- Tanaka did not make the adjustment for the phase (see 
([5|) on Section 2), and in his paper, k is defined as fc = V'o/%/3ai, where is the ratio of the 
incident wave amplitude oq to the reference fiow depth Hq. Tanaka's maximum amplification 
occurs when -00 — 37.8° and — 0.3. Evidently the incident wave angle is too large to make 
the approximation of tanipQ ^ ipQ. Furthermore, the phase adjustment made for the KP limit in 
([5]) cannot be considered negligible. To reflect the finite values of incident angle and to ensure 
the incident wave to be a KdV line-soliton, one must use k = tan ipo j y/2Ao where Aq is defined 
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Fig. 8: Amplification growth a of the stem wave for the conditions of t/jq = 30°, and - -x - -, 
Ao = 0.086; - - o - -, = 0.108; - - - -, = 0.161; = 0.212; = 

0.312; - - A - -, ^0 = 0.413. 

Table 1: Amplification factor a of the stem waves for different Aq {k = tan ■ipo/\/2Ao) with 
tpo = 30°: Q;x=7i.i(Exp.) are the laboratory data a,t X = 71.1 and ax=7i.i(KP) are calculated 
from the corresponding KP exact solutions. ajf=oo(Exp.) are estimated from the exponential 
curve fit to the laboratory data, and q;x=oo(KP) are from (23) and (24). In the last row of 
Aq = 0.413, the values of a in the brackets are obtained a.t X — 50.8, because of the wave 
breaking immediately after this point; hence, the greater amplification cannot be realized. 
Notice the large deviations in the estimates Q;x=oo(Exp.) in the boxes from the theoretical 
predictions for the cases near k = 1. 



Ao 


k 


Q!X=71.l(Exp.) 


Qx=7i.i(KP) ax-- 


=oo(Exp.) ax 


= oo(KP) 


0.086 


1.392 


2.10 


2.36 


2.13 


2.36 


0.108 


1.242 


2.13 


2.51 


2.19 


2.51 


0.161 


1.017 


2.24 


3.38 


2.33 




3.38 




0.212 


0.887 


2.33 


2.43 


2.46 




3.56 




0.312 


0.731 


2.52 


2.61 


2.92 


2.99 


0.413 


0.635 


(2.48) 


(2.54) 


3.94 


2.67 



by (10), that is, the /c-parameter should be 

tan ipQ 



^/3{aQ/ho)cos V'o 



(25) 



As we discussed in Section 2, the wave amplitude Uq in this way represents that of the cor- 
responding incident line-soliton, i.e. = do/ft-o- The difference in the value of k is significant 
when the value of iJjq is not small, obviously. 



The comparison with the numerical results given by Tanaka [27] is shown in Figure 9a 
Note that his numerical experiments were performed for a wide range of the parameter k by 
changing the angle of incidence "00 for a finite amplitude wave with = 0.3. His numerical 
experiments fails to achieve the critical amplification of 4 at fc = 1.0. Also plotted in Figure 



9a are the present laboratory data of which parameters are comparable with Tanaka's data: 



we select the data taken at the farthest location X — 71.1 with the incident wave amplitude 
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fli = 0.277 and incident wave angle ipo = 20°, 30°, and 40° (or equivalently, — 0.367, 0.312, 
and 0.244, respectively). The laboratory data are in good agreement with Tanaka's numerical 
results. It is noted that the observed stem- wave amplitudes are still growing at X = 71.1 as 
shown in Figure [Sj 

Unlike the limited spatial domain constrained by the laboratory wave tank, numerical sim- 
ulations do not have as severe limitations in space and time as the laboratory environment: 
Tanaka's data [27 are almost asymptotically stable a.t X = 150. Although there are slight dif- 
ferences in values of Aq and X, Figure [9a| demonstrates that the present laboratory results are 
in good agreement with Tanaka's numerical results: amplification of the numerical predictions 
are slightly greater than the laboratory data, because of the greater incident wave amplitude 
and the asymptotic amplification in Tanaka's numerical results. Most importantly, Tanaka's nu- 
merical results appear to be in good agreement with ( 23 1 and ( 24 1 of the KP theory as shown in 
Figure [9a| except near the critical point k = 1.0, in spite of the large incident wave angles ipo- It 
is emphasized that our interaction parameter k based on the KP theory (fc = tanipo/y/2Ao) is 
the one that leads to significant improvement in the comparison. As discussed earlier, T ana ka's 
numerical data were significantly deviated to the lower values of k, because instead of (25), he 
used Miles's parameter k = ^o/V^(H defined in [1T],[52]. While ta.n'ipo/y/2Ao w ■0o/v3ai for 
infinitesimal ipQ, it is critical to make proper interpretations for the approximated parameters 
when the mathematical theory is compared with the experimental data with finite values of ipQ. 



None the less, the good agreement shown in Figure 9a is surprising for the large k with large 
incident wave angles. One should emphasize that those results for large k values fit better in 
the KP theory, but not in the theory of regular non-grazing reflection described in [5],[27]- In 
fact, Funakoshi's results [8l with a small amplitude incident wave are in excellent agreement 
with the KP theory for all k values with this new interpretation of the /c-parameter. 





(a) 



(b) 



Fig. 9: Comparison of the present laboratory data of the stem-wave amplification with the 
previous numerical results given by Tanaka [27], and — , theoretical prediction of (231 and (24): 
(a) the laboratory data • at X = 71.1 with the similar conditions (t/jq and Aq) as the numerical 
data aX X — 150 marked by A; (b) the present data are at X = 71 A and the incident wave 
angle Vo = 40°, ■; = 30°, •; = 20°, ♦; = 30° at X ^ +oo, o. The A's are the 
numerical results as in (a). 



Figure 9b shows the stem wave amplification measured in the laboratory at X = 71.1 (our 
farthest measuring location) with ipQ — 20° and 30°, and a,t X = 61.0 with ipQ = 40°; also 
plotted are the KP prediction (23), (24), and Tanaka's numerical results. Note that in the case 
of ipo — 40°, the data at X = 71.1 are slightly contaminated due to the flow contraction owing 
to the limited breadth of the wave tank; hence the data a,t X — 61.0 are used instead. Also 
plotted in the figure are the estimated asymptotic amplifications that are extrapolated from the 
laboratory data as shown in Table [T] The present laboratory data with the incident wave angle 
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(b) 

Fig. 10: Solutions of the KP equation for V^o = 30°: (a) at X = 71.1 (T = 41.05), Aq = 0.112, 
the 0-type (left) and Ao = 0.212, the (3142)-type (right); (b) at X = 208 (T = 120), Aq = 
0.112, the 0-type (left) and Aq = 0.212, the (3142)-type (right). Notice that the stem parts of 
both those solutions are similar at X — 71.1, but later the (3142)-solution has longer stem, the 
Mach stem. 



ifiQ — 40° show excellent agreement with (23 1 of the 0-type solution (fc > 1) as shown in Figure 



|9b[ However, as the incident wave angle decreases, the measured amplification deviates from the 



theoretical prediction. It is surprising to observe in Figure 9b that for the (3142)-type (fc < 1) 



the amplification predicted by ( 24 1 should increase as the value of k increases; our laboratory 
results show opposite. The laboratory data show no sign of transition from the (3142)-type to 
the 0-type reflection. All the present laboratory results appear to exhibit the characteristics of 
Mach reflection, i.e. a) the reflected wave angle is larger than the incident angle; b) the stem 
length continually grows. The smaller the incident wave amplitude, the smaller the amplification 



factor a of the stem wave, and the farther the deviation from the theoretical prediction of ( 23 ) 
for fc < 1. 

A possible explanation for the large discrepancy that appears in the case of V^o = 20° may be 
due to the limitation of propagation distance in the laboratory. With a small incidence angle, 
the development of Mach refiection may take a very long distance to reach its asymptotic 
state. Examination of our KP solutions supports this. Figure |10a| shows the KP solutions at 
X = 71.1 based on the recipe provided in Section 2 for the cases with Aq = 0.112 and 0.212 
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with V'o = 30°. Note that the former case Aq = 0.112 is the 0-type because ipo > tpc = 25.3°, 
and the latter case = 0.212 is the (3142)-type because ipQ < ipc = 33.1°. The same waves but 
at X = 208 are presented in Figure [TOb) At X = 71.1, the stem length of the 0-type solution 
appears similar to that of the (3142)-type, while they are distinctively different a,t X — 208. 
This observation indicates that both the 0-type and the (3142)-type would grow their stems 
at the similar rate at least up to X = 71.1. However, the growth of the stem is ceased for the 
0-type once it reaches its equilibrium state. Evidently, this is not the Mach stem formation 
described by Miles |22j . On the other hand, the Mach stem of the (3142)-type appears continual 



growth (see Figure 10 right). Most of our laboratory data were recorded at X = 71.1, where 



both 0-type and (3142)-type are still being developed. The KP solutions in Figure 10 suggest 
that the laboratory tank must be too short to distinguish these two types of the KP solutions. 

Table [2] shows the KP predictions of amplification a with tpQ = 20° at X = 71.1; the 
predicted amplifications are in excellent agreement with the measured values. We extended 
the KP predictions to a longer distance, X = 350, close to the asymptotic amplifications. 
The excellent agreement in amplification of the predictions with the measurements at X = 
71 A suggests that the shortcoming of the laboratory experiments should be attributed to the 
available propagation distance that must be insufficient to establish its asymptotic condition. 



Table 2: Asymptotic amplification of the stem waves with ipo = 20°: k = tan tpo/^/2Ao, and 
Q!j)(:=7i.i(Exp.) is the measured stem-wave amplification. ax=7i.i(KP) is the KP prediction at X 
= 71 A, and ax=35o{^P) is that aX X = 350 where the reflection should become its asymptotic 
state. 



^0 


k 


ax=7i.i(Exp.) 


QX=71.l(KP) 


ax=35o(KP) 


0.127 


0.722 


1.89 


1.84 


2.87 


0.189 


0.591 


1.95 


1.93 


2.52 


0.249 


0.516 


1.99 


2.08 


2.33 


0.367 


0.425 


2.01 


1.99 


2.03 



5 Summary and Conclusions 

Amplification of obliquely incident solitary wave onto a vertical wall was analyzed based on 
the KP equation. The KP equation is equivalent to the approximation of the theory developed 
by Miles Consequently the same predictions result for the asymptotic conditions, 

which include the maximum four-fold amplification at the critical state. Unlike Miles's theory, 
however, the KP equation enables us to analyze the development processes. Our theoretical 
predictions were analyzed with our laboratory experiments as well as the previous numerical 
experiments by Tanaka (27]. 

Careful attention must be taken when theoretical predictions are compared with laboratory 
or numerical data. Adequate adjustments must be made to link the physical quantities with 
the mathematical parameters. Specifically, the interaction parameter k should be computed 



with (25), i.e. k = taniljo/V^A^ — tan ijjo/{y^5ai cos V'o)j instead of fc = ipo/y/Soi. The former 
parameter k ensures that the incident wave is a KdV soliton, while the later would not. It is 
because the incident wave angle ipo is not infinitesimal but finite in the experiments. Discrepancy 
in amplification reported i n T anaka' s paper [27' improves significantly with the use of proper k 



as demonstrated in Figure 9a We also confirm that Funakoshi's results 8J with smaller incident 
amplitudes are in an excellent agreement with the KP prediction with our interaction parameter 
k. 

The present laboratory results of wave amplification are consistent with Tanaka's numerical 
results if and only if the data with similar conditions are compared. However, when the compar- 
ison is extended to include all other data, the present results turn out different from both KP 
theory and Tanaka's numerical results. The KP theory indicates that our laboratory apparatus 
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is too short to demonstrate the difference between the 0-type and the (3142)-type solutions. 
In fact, the stem wave along the wall grows at the similar rate for both 0-type and (3142)-type 
along the available length in the laboratory. The theory predicts that the growth rate soon 
deviates beyond the distance of the apparatus the 0-type solution ceases the growth when it 
reaches its equilibrium state, while the stem wave continually grows for the (3142)-type as the 
characteristic of the Mach stem. Considering a long distance (time) for an obliquely incident 
soliton to become the asymptotic state, especially for a small incidence tpQ, it would be difficult 
to realize the theoretical predictions in a real-fluid laboratory environment. 

This work is supported by Oregon Sea Grant Program No. NA06OAR4170010. One of the authors 
(Y.K.) is partially supported by NSF Grant DMS-0806219. 
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